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‘‘Disentangling nestedness’’ disentangled 


ARISING FROM A. James, J. W. Pitchford & M. J. Plank Nature 487, 227-230 (2012) 


Analytical research indicates that the ‘nestedness’ of mutualistic net- 
works facilitates the coexistence of species by minimizing the costs of 
competition relative to the benefits of facilitation’. In contrast, James 
et al? recently argued that a more parsimonious explanation exists: 
the persistence of a community and its constituent species depends 
more on their having many interactions (high connectance and high 
degree, respectively) than for these interactions to be organized in any 
particular manner. Here we demonstrate that these conclusions are an 
unintended consequence of the fact that the methodology of ref. 2 
directly changed the number of interactions of each species—and hence 
their expected persistence. When these changes are taken into account, 
we find a significant, positive relationship between nestedness and 
network persistence that reconfirms the importance of nestedness in 
mutualistic communities’’. There is a Reply to this Brief Communi- 
cation Arising by James, A., Pitchford, J. W. & Plank, M. J. Nature 500, 
http://dx.doi.org/10.1038/nature12381 (2013). 

Given a network, one can robustly quantify the relative numbers of 
specialist to generalist species via the degree distribution**. A network’s 
degree distribution is of considerable importance, because studies have 
repeatedly highlighted the significant, positive relationship between a 
species’ number of mutualistic partners and its survival probability’ **. 
This distribution alone is also capable of driving many higher-order 
network properties’, not to mention the fact that the degrees of species 
are phylogenetically constrained themselves*. For these and other rea- 
sons, studies across the ecological-network literature**” have empha- 
sized the need to take the degree distribution into consideration when 
assessing the significance of the myriad patterns observed in nature”""’. 

Unfortunately, when comparing empirically observed networks to 
random networks, the authors of ref. 2 seem to have overlooked this 
critical link between changes in the degree distribution and species’ 
survival. As a direct consequence, the specialists in their random net- 
works became less specialist and the generalists less generalist’. Yes, 
the random networks were observed to be more persistent (Fig. 1a), 
but this was not in fact an indication that nestedness is unimportant’. 
Instead, this increase in persistence was a result of the random net- 
works having more homogeneous degree distributions”'’, and that 
the most vulnerable species in the empirical networks almost always 
had more interactions in the corresponding randomizations. Here 
this distinction is of critical importance because species’ degrees are, 
in fact, “a better predictor of individual species survival”. “The more 
the merrier” indeed”. 

To quantitatively validate these results, we repeated a key analysis 
of ref. 2 to measure the relationship between nestedness and persist- 
ence while paying explicit attention to changes in the network’s degree 
distribution (Methods). On taking the small but critical step of con- 
trolling for the increased homogeneity of the degree distributions, we 
observe a significant, positive relationship between nestedness and 
persistence (Fig. 1b). In addition, we reach the same conclusion whether 
we account for changes in the degree distribution statistically or by 
repeating the analysis while generating the randomized networks with 
a null model that explicitly maintains the observed degree distribution 
(Fig. 1c, Methods and Appendix). All else being equal, our results here 
illustrate that, the greater the nestedness of a community, the greater 
indeed is that community’s persistence. 

Given an observed number of species and interactions in a com- 
munity, a prevailing question across the ecological literature is whether 
or not some ways to structure those interactions (for example, nested- 
ness) lead to more persistent communities. Although the number of 
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Figure 1 | Within our regression analysis, the relationship between 
nestedness and persistence in mutualistic networks depends integrally on 
changes in the degree distributions of the networks. a, If these distributions 
are allowed to change but are uncontrolled for, nestedness appears to be 
negatively correlated to persistence (P < 10“). b, c, However, when these 
changes are appropriately controlled for—either statistically (b) or in the null 
model for randomization (c)—there is a significant positive relationship 
between nestedness and persistence (P< 10 *andP<10 “4, respectively). The 
same general conclusions reached here for the probabilistic null model hold for 
other, non-degree-preserving randomizations’. 


mutualistic interactions of a species plays an important role in its 
survival***"*, we find unambiguous support for the added importance 
of the way in which mutualistic interactions are organized—the true 
architecture of biodiversity'*. Echoing ref. 2, our findings re-emphasize 
the importance of carefully considering the interplay between all 
potential sources of variation™ in ecological models. Otherwise, one 
runs the risk of further entangling models that are sufficiently tangled 
already. 


Methods 


For 59 empirical networks, we generated 250 randomized networks and for each we 
simulated persistence (the fraction P of surviving species in each simulation) across 
250 parameterizations of a dynamic mutualistic model’. We quantified the relation- 
ship between persistence and nestedness with a mixed-effects logistic regression’ 
that takes the form logit(Pj,) = Bo + BiM; + BoC; + P3 Wij + BaNy + ni + rij + Eijk 
Here the indices i, j and k indicate the empirical network, network randomization 
and model parameterization, respectively, fo is a constant, the slopes f1, P2, B3 and 
ba quantify the importance of network magnitude’ M, connectance’ C, relative 
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degree homogeneity’? W, and nestedness” N, respectively, the random effects n; 


and rj control for variance across networks and randomizations, and ¢; is the 


model residual. Variance inflation factors gave no indication of multicollinearity in 
this model. 
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Appendix 

We randomized the empirical networks with two null models: the probabilistic and 
fixed (or swap) algorithms’. For our purposes here, the key distinction between the 
two is that the probabilistic model generates random networks with quantitatively 
more homogeneous degree distributions than those observed empirically (Wj > 0) 
whereas the degree distribution is strictly conserved in networks generated by the 
fixed model (Wj = 0). The statistical analyses presented here were performed in R 
version 2.15.3 (http://R-project.org/) using the glmer function in package Ime4 
version 0.999999-0 (http://Ime4.r-forge.R-project.org). Code to perform the network 
randomizations and dynamic simulations in Matlab (http://www.matlab.com/) 
and the mixed-effects logistic regressions in R (http://R-project.org/) is available 
from the Dryad Digital Repository at http://dx.doi.org/10.5061/dryad.p2gq8. 


REPLYING TO S. Saavedra & D. B. Stouffer Nature 500, http://dx.doi.org/10.1038/nature12380 (2013) 


Saavedra and Stouffer’ claim that the results of James et al.? are a 
consequence of the method used to randomize interaction matrices. 
We recognize the importance of examining alternative randomization 
schemes and have repeated our analysis using their methods. 
However, we find no evidence that ‘reconfirms the importance of 
nestedness in mutualistic communities”. 

Repeating the analysis of figure 2 in ref. 2, using the swap rando- 
mization scheme’, which does not change degree distribution, con- 
firms our finding that the persistence of real networks is not related to 
their nestedness. Although more of the empirical networks are less 
nested than their randomized counterparts under this scheme, con- 
trary to the accepted result of ref. 4, there is no useful correlation 
between nestedness and persistence (Fig. 1). Therefore, the ‘small 
but critical step’ of accounting for degree heterogeneity does not 
produce a positive relationship between nestedness and persistence. 

The results in figure 1 of ref. 1 represent relationships between 
nestedness and persistence among randomizations of individual net- 
works. They do not imply that, given two observed networks, the 
more nested network is more likely to have the higher persistence 
as claimed in ref. 5. We have performed the general linear mixed 
model (GLMM) analysis advocated in ref. 1. This shows that >90% 
of the variance comes from variance between groups (networks) and 
<10% comes from variance within groups. This highlights the lack of 
consistency across the groups, and that any effect of nestedness is 
dwarfed by the random effects of the GLMM. 

The NODF* definition of nestedness used in refs 1 and 2 is one of 
several possible metrics. For example, the nestedness metric used in 
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Figure 1 | Accounting for degree distribution does not give a meaningful 
relationship between nestedness and persistence. Repeating the 

results of figure 2 in ref. 2 using the swap randomization scheme, 

which does not change connectance or degree distribution, there is no 
useful correlation between the change in nestedness (relative to the 
empirical network) and the change in the persistence of the dynamic 
model. Each point represents the average of 100 randomizations of an 
empirical network. 
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ref. 5 is invariant under the swap randomization scheme of ref. 1, 
precluding a GLMM approach. Under the metric of ref. 7, the 
GLMM reveals a negative relationship between nestedness and per- 
sistence. That the conclusions of ref. 1 are sensitive to the choice of 
metric indicates that they cannot be used to draw general conclusions 
about the effects of nestedness. In contrast, the results in ref. 2 are 
robust to the choice of metric. 

Is nestedness important for predicting persistence in these models? 
Our results, confirmed by the methods of ref. 1, show that it is less 
important than: network size; connectance; degree distribution; 
intrinsic growth rates; competition coefficients; and the strength of 
the mutualistic interactions. If two ecosystems can be found that share 
all these properties then, under the specific dynamic model tested 
here, the more nested ecosystem may (depending how nestedness is 
defined) be more likely to have a higher persistence. However, if any of 
these properties differ between the two ecosystems, then any effect of 
nestedness is likely to be unimportant. 

In conclusion, nestedness is an interesting abstract network prop- 
erty that undoubtedly influences the statistical behaviour of large 
systems of differential equations’. However, general conclusions 
allowing nestedness to be used as a predictor of empirical biodiversity 
cannot currently be justified. 
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